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ABSTRACT 

We analyse N-body galaxy merger experiments involving disc galaxies. Mergers of disc- 
bulge-halo models are compared to those of bulge-less, disc-halo models to quantify 
the effects of the central bulge on the merger dynamics and on the structure of the 
remnant. Our models explore galaxy mass ratios 1:1 through 3:1, and use higher bulge 
mass fraction than previous studies. A full comparison of structural and dynamical 
properties to observations is carried out. The presence of central bulges results in 
longer tidal tails, oblate final intrinsic shapes, surface brightness profiles with higher 
Sersic index, steeper rotation curves, and oblate-rotator internal dynamics. Mergers of 
bulge-less galaxies do not generate long-lasting tidal tails, and their strong triaxiality 
seems inconsistent with observations; these remnants show shells, which we do not 
find in models including central bulges. Giant ellipticals with boxy isophotes and 
anisotropic dynamics cannot be produced by the mergers modeled here; they could be 
the result of mergers between lower luminosity ellipticals, themselves plausibly formed 
in disc-disc mergers. 

Key words: galaxies:interactions- kinematics and dynamics- structure- elliptical - 
numerical simulation 



1 INTRODUCTION 



Hierarchical galaxy formation models assume a merger ori- 
gin for spheroids (e.g. Cole et al. 2000). Such merger hy- 
pothesis was initially put forward by Toomre (1977), and 
Toomre's formulation is still useful for detailed studies of 
the links between merger dynamics and elliptical galaxy 
structure. Observations of merger sequences (e.g. Hibbard 
& Mihos 1995) generally support the merger hypothesis, as 
mergers can explain most of the characteristics of elliptical 
galaxies deemed 'peculiar', such as shells, ripples, kinemati- 
cally distinct cores, misaligned rotation axes and disky-boxy 
isophotal shapes (see Schweizer 1998 for a review). N-body 
simulation techniques are ideal to check whether the col- 
lisionless evolution of galaxies in mergers leads to objects 
showing the structural and kinematic properties of ellipti- 
cals, and many papers have studied the merger hypothesis 
using these techniques. An underlying question, which re- 
mains open today, is whether the solution to the divergences 
between models and real ellipticals requires the refinement 
of the N-body models; the modeling of new stars formed out 
of the gas component of the precursor galaxies; or the ruling 
out of the merger hypothesis in favor of classical collapse 
models. 

Since the late eighties, improvements in the models and 
the computing power have led to a progressive increase in the 
understanding of the merger process. Barnes (1988, 1992) 



used equal-mass model galaxies composed of a non-rotating 
bulge, a disc and a dark matter halo. The remnants from 
his experiments show de Vaucoulers profiles compatible with 
real E's and they are supported by random motion rather 
than by rotation. Hernquist and collaborators (1992, 1993; 
Hernquist et al. 1993; Heyl et al. 1994, 1995, 1996) studied 
equal-mass mergers of both bulge-less systems and systems 
with a bulge. Bulge-less systems were deemed to be unlikely 
progenitors of today's elliptical galaxies due to the low den- 
sity in the inner regions of discs. This difficulty was first 
put forward by Ostriker (1980) in a short but famous paper, 
and later by Carlberg (1986). Two ingredients may in prin- 
ciple alleviate this problem: the introduction of gas, or of 
a third component, namely a bulge. Hernquist et al. found 
that, as expected, the presence of bulges in the progenitors 
boosts the central density of the remnant, thus providing a 
plausible solution to the central density problem. 

Barnes (1999) extended the exploration of N-body 
mergers to unequal-mass progenitors. His mergers with 3:1 
mass ratio generally lead to fiattened systems, irrespective 
of the orbital configuration, in contrast with 1:1 merger rem- 
nants, whose shapes strongly depend on the initial orbital 
geometry. This led to the proposal that SO galaxies may be 
the result of mergers between unequal-mass discs. 

The link between unequal-mass mergers and remnants 
harboring discs was furthered by Naab, Burkert & Hernquist 
(1999), who analysed the isophotal properties of 1:1 and 3:1 
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mass-ratio mergers. Isophotes from 1:1 remnants tend to 
have boxy deviations from elliptical shape, while 3:1 merg- 
ers lead to disky objects instead. They thus showed that 
coUisionless mergers can account for the range of isophotal 
shapes of elliptical galaxies, without recourse to gas dynam- 
ics to explain disky objects. More recently Naab & Burk- 
ert in a series of papers (2001, 2003) extended the study 
with a large number of disc-disc mergers similar to the ones 
found in previous papers. They find that mass ratio and ini- 
tial orbital configurations are important in the shaping of a 
number of characteristics of the remnants. 

To what degree kinematic diagnostics support the 
merger hypothesis is unclear at present. Reproducing the 
kinematic peculiarities in ellipticals, such as misaligned ro- 
tation and counter-rotation, with coUisionless mergers is 
not difficult; we (Balcells & Gonzalez 1998) modeled col- 
lisions with mass ratios 1:1, 2:1 and 3:1, using galaxy mod- 
els which included bulge components in the progenitors. 
Unequal-mass mergers in retrograde orbits lead to remnants 
with counterrotating cores, without recourse to dissipative 
gas dynamics (Hernquist & Barnes 1991) to concentrate the 
counterspinning material in the nuclear regions of the rem- 
nant; also, Bendo & Barnes (2000), using merger simulations 
with mass ratios 1:1 and 3:1, found that the merger models 
can explain kinematic properties of real ellipticals such as 
misaligned rotation and counter-rotation. However, whether 
such mergers can explain the more fundamental question of 
the amount of rotation of ellipticals remains to be answered. 
Part of the difficulty rests in whether the goal is to reproduce 
the slow rotation of giant ellipticals or the rapid rotation 
of intermediate- and low-luminosity ellipticals, and whether 
the comparison is to be performed in the inner parts of the 
galaxies or beyond the effective radius, where observational 
data is scarce but with some evidence for rapid rotation 
(Rix, CaroUo, & Freeman 1999). Bendo & Barnes (2000) 
claim that their models reproduce the rapid rotation found 
in the sample of Rix et al. (1999). However, Cretton et al. 
(2001), using similar simulations, conclude that their models 
are unable to yield the rapid rotation in the Rix et al. data. 

The study of the asymmetry term /13 in the line-of-sight 
velocity distributions (LOSVD) provides additional clues. 
Observationally, besides the well-known trend that hs has 
opposite sign to the mean velocity (e.g. Bender et al. 1994), 
elliptical galaxies show high hs values near the galaxy cen- 
tres, where v/a is low (Bender et al. 1994; their Figure 15). 
Naab & Burkert (2001) tried to reproduce this result using 
mergers with 1:1 to 3:1 mass ratios. They concluded that 
their merger remnants do not follow the trends unless they 
added a disc with at least 10 to 20% of the spheroid mass,. 
They proposed gas as the means to produce this disc. 

Thus, although mergers of disc systems account for 
a great deal of both global and fine-structure features of 
ellipticals, several issues remain open. We highlight four 
questions: (i) can the central densities of ellipticals be ex- 
plained through coUisionless merger dynamics of the precur- 
sors' bulges, or are they the traces of gas infall with ensuing 
star formation, (ii) do the diskiness properties of ellipticals 
require gas, (iii) what conclusions can be drawn from the rel- 
atively unexplored 3:1 mergers, and (iv) do the rotation and 
LOSVD properties of coUisionless merger remnants match 
those of real galaxies. 

In the present paper we address these questions through 



Table 1. Initial input parameters to generate dbh and dh models 
with the GalactlCS code. 



model 


# halo 










dbh 


30000 


-4.6 


1.42 






dh 


30000 


-3.0 


1.5 






model 


#disc 


Md 


Rd 


-^^outcr 




dbh 


12000 


0.867 


1 


5 


0.1 


dh 


15000 


1.5 


1 


5 


0.1 


model 


# bulge 


Ph 


o-b 






dbh 


6000 


14.45 


0.714 






dh 
















the analysis of 18 simulations of mergers between disc gala^x- 
ies. We run each merger with initial models comprising disc, 

bulge and dark halo, as well as bulge-less models compris- 
ing a disc and a halo only. We have explored three differ- 
ent mass ratios of the progenitor systems (1:1, 2:1, and 3:1) 
and three different orbital orientations, keeping the initial 
orbits for the 18 runs the same. The survey is clearly not 
comprehensive, but it expands on previous works in sev- 
eral ways. The studies by Barnes (1992), Hernquist (1992), 
Hernquist (1993), Hernquist et al. (1993a), Hernquist et al. 
(1993b), Heyl et al. (1994), Heyl et al. (1995), Heyl et al. 
(1996) involve equal-mass galaxies, hence do not take into 
account the effects of different mass ratios on the struc- 
ture of the final system. Barnes (1999), Balcells & Gonzalez 
(1998), Naab et al. (1999), Bendo & Barnes (2000), Cretton 
et al. (2001) and Naab & Burkert (2001, 2003) do explore 
mergers between progenitors with mass ratios of 1:1 and 
3:1 including a bulge. Naab & Burkert (2003) extend their 
range up to mass ratio 4:1. Our models are complementary 
to those above in that our galaxy models have more mas- 
sive bulges, with disc-to- bulge ratio D/B=2:l, i.e. we model 
'earlier type' galaxy mergers. Finally, we analyse bulge-less, 
non-equal mass mergers, which are not addressed in previ- 
ous papers. Together, our 'early-type' and bulge-less models 
provide a wide base-line in bulge-to-disc ratio (B/D), which 
is useful to address the central density problem as well as 
the controversy on the rotation properties of the merger rem- 
nants. 

The initial models and the setup of the simulations 
are described in Section 2. Section 3 presents a comprehen- 
sive analysis of the merger remnants, including tidal tails, 
shells, intrinsic shapes, surface density profiles, isodensity 
contour shapes, rotation curves, pressure support, and kine- 
matic anisotropy. In section 4 we study the effect of particle 
resolution. We discuss our results in Section 5. 



2 INITIAL MODELS 

The initial models were built using the Kuijken-Dubinski 
GalactlCS code (Kuijken & Dubinski 1995), made freely 
available by the authors. The halo is built from an Evans 
model (Kuijken & Dubinski 1994); the bulge is a King (1966) 
model, thus having a nearly-exponential surface density pro- 
file; and the disc is a Shu (1969) model generalized to three 
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Table 2. Initial mass ratios and disc central radial velocity dis- 
persion for dbh and dh models. 



model 


(Afluni/A^tot)Bo„t 


(Mum/Mot )total 


O"R.0 


dbh 


0.459 


0.241 


0.410 


dh 


0.463 


0.129 


0.320 



dimensions; its surface density is nearly exponential, and it 
follows a sech^(z) profile vertically. Parameters for the ini- 
tial realizations are given in Table Haloes are given by 
the global central potential V'o and the asymptotic circular 
velocity vq. Discs are given by their mass Ma, the exponen- 
tial scale length R^, the outer truncation radius Routa and 
the vertical scale height z^- Bulges are given by their central 
density pb and central velocity dispersion ah- The number 
of particles and the mass of each component are given in 
columns (2) and (3), respectively (units are described in 
section im . Neither the bulges or the haloes are given net 
rotation. It should be noted that the models generated by 
the code where scaled to the desired masses as presented 
below. Table|5|lists the ratios of luminous-to-dark masses in 
each of the models, and the central radial velocity dispersion 
for the models with a Mtot = 1- 

Our disc-bulge-halo model (hereafter dbh models) , iden- 
tical to model A from Kuijken & Dubinski's (1995) models 
for the Milky Way, has a bulge-to-disc mass ratio B/D=l:2. 
Basic properties of this model are shown in Figure left 
panels. The circular velocity curve is flat out to the disc 
truncation radius (iJoutcr = 5) and gently declines further 
out. Note that the model does not have a maximum disc, 
rather it is bulge-dominated out to r = 1.5f?d and halo- 
dominated elsewhere. The surface density profile shows an 
exponential profile for the bulge component, and a charac- 
teristic exponential profile in the disc region. Note however 
that the disc is not exponential all the way to the galaxy 
nucleus. 

The bulge-less, disc-halo models (hereafter dh), were 
also built using the Kuijken & Dubinski (1995) algorithm. A 
bulge-less disc is quite unstable to bar formation (Hernquist 
1992), and building a stable model proved to be difficult. We 
obtained a solution by lowering the concentration of the halo 
and rising that of the disc. The model has nearly maximum 
disc, a low Toomre parameter Q=1.4 at r ~ 3.5i?d, and a 
surface density profile that deviates from a pure exponential 
in the inner regions (Fig.0 right panels). 



2.1 Units 

We use model units so that the constant of gravity is G = 1. 
The total mass of the smallest system is taken to be M = 1. 
The scale length of the smallest disc is Rd — 0.2. A set of 
physical units that match the dbh models to the Milky Way 
are: 




0.6 
0.5 
0.4 
S 0.3 
S 0.2 
0.1 
0.0 
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Figure 1. Left panels: the initial dbh models. Right panels: 
the initial dh models. Top panels: Circular velocity curves (solid 
lines), and contributions from the disc (dashed), the bulge (dot- 
dashed), and halo (dotted). Middle panels: Toomre Q parame- 
ter. Bottom panels: face-on surface density profiles for the lumi- 
nous matter (solid lines), for the disc (dotted) and for the bulge 
(dashed). is the exponential scale length of the disc. 



[v] = 315 km/s. (4) 



2.2 Stability of initial models 

The stability properties of the initial dbh and dh models 
are analysed by evolving the model galaxies in isolation and 
studying the time evolution of the radial mass distribution, 
the disc thickness and the circular velocity curve. These runs 
last for about 25 half-light crossing times of the disc galaxy. 

The results for dbh models are summarized in Figure |5| 
The radial mass distribution (Fig.|2Ji.) and the rotation curve 
Vc = \/ GM(r)/r (Fig.|5|D) remain stable. Figure shows 
the final circular velocity curve, to be compared with the 
one in Figure Q The disc vertical scale height increases with 
time due to bombardment by the halo particles (Fig.|5J;). 
The increase of ~70% over the duration of the experiment is 
comparable to that measured by Kuijken & Dubinski (1995). 

Figure El shows the stability results for model dh. Over 
35 half-light crossing times the system suffers a mild reorga- 
nization of the mass distribution. The disc shows some tran- 
sient spiral features but is stable against bar formation. Fig- 
ure|2Ji shows the evolution with time of fractional mass radii. 
We find that after relaxation the curves have not changed 
much. Figure shows that the disc thickness increases by 
60%. 



[M] = 3.24 X lO" Mq, 

[L] = 14.0 kpc, 

[T] = 4.71 X lO'^ yr, 

with 



(1) 
(2) 
(3) 



2.3 Scaling rules 

We scale the sizes of model galaxies of different masses using 
the TuUy-Fisher (TF) relation, assuming constant mass-to- 
light ratios as detailed below. According to the TF relation. 
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Figure 2. (a) Evolution of the radial luminous mass distribution 

for the initial dbli model, (b) The final circular velocity for the 
total mass (diamonds), the disc (squares), bulge (triangles) and 
halo (crosses). The stability can be check comparing with the 
initial curves for the total mass (solid line) , the disc (dashed line) 
the bulge (dot-dashed line) and for the halo (dotted line), (c) 
Evolution of the thickness of the disc during the stability run. 



Figure 3. (a) Evolution of the radial luminous mass distribution 

for the initial dh model, (b) The final circular velocity for the 
total mass (diamonds) the disc (squares) and the halo (crosses). 
This can be compared with the curves for the initial model for the 
total mass (solid line), for the disc (dashed line) and for the halo 
(dotted line) . (c) Evolution of the thickness of the disc during the 
stability run. 



the maximum of the velocity curve scales as a power a of 
the luminosity of the galaxy: 



L OC (Knax)", 



(5) 



where a ranges from 2.5 (Ziegler et al. 2002) for the B-band 
for distant galaxies, to 3.2 for the B band, 3.5 for R and 4 

for the I band (Sakai ct al. 2000) for nearby galaxies. 

Now, assuming virial equilibrium, 2Et: = —Ew, we 

have 



V OC 



(6) 



where is twice the kinetic energy per unit mass of the 
galaxy and re is the gravitational radius, defined by the 
equation Ew ~ —GMtot/rc- 

Assuming that the luminous mass scales with the lumi- 
nosity {Mium/L ^constant), 



Mum tX (Vmax)", 

with V (X Vmax we obtain: 



(7) 



(8) 



Finally, we assume that Mtotai/A/ium =constant, and 
we obtain the scaling law: 
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Figure 4. Luminosity versus maximum rotation velocity for the 
initial disc models. The luminosity is taken equal to the lumi- 
nous mass (bulge+disc in dbh models and disc in dh models). 
Diamonds are dbh models. Solid line is a fit to those points. The 
slope of this line is 3.51. Squares are dh models and the dotted 
line is a fit to these points; it has a slope of 3.45. 



We have chosen a — 3.5, matching the TF exponent in 
the R-band (Sakai et al. 2000). 

Figure|l| shows that the models scaled with this pro- 
cedure indeed follow a Tully-Fisher relation. The observed 
shift between the two relations reflects a different zero-point 
due to the different masses of the luminous components of 
dh and dbh models. 

2.4 Merger experiments 

We run mergers involving either two dbh models or two dh 
models. Mass ratios of 1:1, 2:1 and 3:1 are done in each case. 
For each combination of initial model and mass ratios, three 
relative orientations of the galaxy spins with the orbital an- 
gular momentum are run, yielding a total of 18 independent 
merger models. Merger orbital configurations are not varied 
for each combination of initial model and mass ratio. Hence, 
our experiments are suited for an investigation of the effects 
of mass ratio, of spin orientation and of the presence of a 
bulge on the structure of the final remnant. However, we do 
not intend to map the complete parameter space. 

Merger parameters are listed in Table|3| A right-handed 
coordinate system has been adopted in which the initial 
models are placed on the x-axis, positive toward the less 
massive model and the orbital angular momentum defines 
the positive z-axis. Column (1) gives the model name. The 
number denotes the mass ratio, as given in column (2); note 
that masses reflect total, luminous plus dark components, 
that M2 denotes the most massive system, and that Mi = 1 
in all cases. Model names without an asterisk denote dbh 
mergers, names with an asterisk denote dh mergers. The 
letter describes the spin orientations: 'p' models have a pro- 
grade coupling of both galajcy spins with the orbit; 'a' ("an- 
tiparallel") models have a prograde coupling of the spin of 
the most massive galaxy with the orbit, and a retrograde 
coupling of the least massive galaxy; and V ("retrograde") 
models have both spins retrograde with respect to the orbit. 



Table 3. Input configurations for disc- bulge- halo models (first 9 
entries) and for disc-halo models (maked with *). 



Mod. §i (ei,«!.i),(e2,<^2) n=a e rp„i Tp^, 
(1) (2) (3) (4) (5) (6) (7) 



Ip 1:1 (10,-10);(70,30) 11.25 0.7 3.38 23.23 
la 1:1 (130,60);(70,30) 11.25 0.7 3.38 23.23 
Ir 1:1 (170,-10);(110,30) 11.25 0.7 3.38 23.23 



2p 2:1 (10,-10);(70,30) 24.76 0.7 7.43 61.94 

2a 2:1 (130,60);(70,30) 24.76 0.7 7.43 61.94 

2r 2:1 (170,-10);(110,30) 24.76 0.7 7.43 61.94 

3p 3:1 (10,-10);(70,30) 26.98 0.7 8.07 60.99 

3a 3:1 (130,60);(70,30) 26.98 0.7 8.07 60.99 

3r 3:1 (170,-10);(110,30) 26.98 0.7 8.07 60.99 



Ip* 1:1 (10,-10);(70,30) 12 0.7 3.60 25.60 

la* 1:1 (130,60);(70,30) 12 0.7 3.60 25.60 

Ir* 1:1 (170,-10);(110,30) 12 0.7 3.60 25.60 



2p* 2:1 (10,-10);(70,30) 20 0.7 6.00 44.97 

2a* 2:1 (130,60);(70,30) 20 0.7 6.00 44.97 

2r* 2:1 (170,-10);(110,30) 20 0.7 6.00 44.97 

3p* 3:1 (10,-10);(70,30) 22 0.7 6.60 44.93 

3a* 3:1 (130,60);(70,30) 22 0.7 6.60 44.93 

3r* 3:1 (170,-10);(110,30) 22 0.7 6.60 44.93 



The exact orientations of the initial spins are given in column 
(3), with 9 and (p being the usual polar and azimuthal spher- 
ical coordinates. We avoided exact parallel or anti-parallel 
spin orientations since these are not representative of a gen- 
eral encounter between galaxies. Columns (4)-(7) list the 
initial separation, ellipticity, pericentre separation and time 
to pericentre for the Keplerian orbit with the same initial pa- 
rameters, respectively. The orbits are slightly sub-parabolic, 
and ensure that the mass distributions overlap at pericentre, 
leading the systems to a fast orbital decay. 

After the merger is complete, systems are let to relax for 
at least 40 half-light radius crossing times before analysis. 

2.5 Merger simulation details 

Merger simulations are run using Hernquist's version 3 of 
the TREECODE (Hernquist 1990) on an Ultra Sparc Sta- 
tion. Softening is set to one-fifth of the bulge half-mass ra- 
dius for the smallest galaxy in each simulation, e = 0.02 
for dbh models, and to one-fifth of the disc half-mass radius 
for dh models, e = 0.09. The tolerance parameter is set to 
0.8, and quadrupole terms are included in the force calcu- 
lation. Individual time step calculation is included. Energy 
conservation errors are below 1%. 

We have modeled the bulge with 6000 particles, the disc 
with 12000, and the halo with 30000 particles for dbh mod- 
els. For dh models the disc has 15000 particles and the halo 
has 30000 particles. In the choice of particle numbers we 
took care that the particle masses for the various compo- 
nents were not too dissimilar to minimize particle-particle 
heating. The maximum ratio of particle masses is 7.5 for dbh 
models and 11.6 for dh models. The total particle numbers 
are rather low but they are comparable to Barnes (1992) 
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and Hernquist (1992, 1993) and not so different from later 
work by these authors and others (Heyl et al. 1996, Naab et 
al. 1999). But our choice of particle number is clearly lower 
than more recent work (e.g. Naab & Burkert 2003). As a 
check of the discreteness and resolution consequences of a 
low particle number, we have replicated one of our experi- 
ments (model 2a), increasing the number of particles by a 
factor of 5. This experiment, named 2aG, was run with the 
GADGET code (Springel et al 2001) on a Beowulf cluster. 
The properties of this model are represented with a filled di- 
amond in the relevant figures throughout the paper. We dis- 
cuss it in Section^] We find that model 2aG shows only mi- 
nor differences with respect to its low-A^ counterpart, which 
suggests that the diagnostics reported in the paper for the 
merger models do not suffer from important discreteness or 
resolution effects. 



3 RESULTS 

The 18 merger remnants of our experiments display a wide 
variety of morphological, kinematic and 'photometric' char- 
acteristics. These must reffect the effects of varying mass ra- 
tio, of the spin orientations and of the presence or absence 
of a bulge. 

Figures|S]and|n]show the time evolution of the luminous 
matter in two representative runs. Figure|K| gives the evolu- 
tion of run Ip. Two disc- bulge- halo models collide with their 
spins partially aligned with the orbital angular momentum. 
Prominent tails are created after the first pass through the 
pericentre. These tails are still prominent well after the in- 
ner parts are merged. FigureJS] shows the time evolution of 
model Ip*. Here two dh models merge from the same or- 
bit as in Figure|S] The tails formed after the first passage 
through pericentre are less prominent, and, by the time the 
merger has completed, the tails have disappeared. 

3.1 Tails 

Tail morphology in the nine dbh merger remnants is depicted 
in Figure |7| which shows the large-scale distributions of the 
luminous matter roughly 1 Gyr after the mergers are com- 
plete. The size of each system has been scaled by the initial 
scale length of the largest disc, Rd2 (see eqn.|5J, taken here 
as a reference yardstick for comparing models with varying 
masses and sizes. 

Our models confirm previous work that prograde equal- 
mass mergers yield two tails, one coming from each of the 
precursor galaxies (Toomre & Toomre 1972, hereafter TT72; 
Barnes 1992). However, model Ip shows that tails from pro- 
grade, equal-mass mergers can be quite asymmetric, when, 
as is the case in our models, spins are not exactly aligned 
with the orbit. Also, tails become highly asymmetric, ap- 
proaching one-tail configurations as soon as the mass ratio 
departs from unity (models 2p, 3p). Hence, galaxies exhibit- 
ing one single, prominent tidal tail are not necessarily the 
result of a minor merger. 

All dbh remnants, from both prograde and retrograde 
merger orbits, show tails. From spin-orbit coupling (TT72; 
Barnes 1992; Hernquist 1993), tails were expected mp mod- 
els (both spins aligned with Jorb) and in a models (primary 
aligned with Jorb), those of the p model being longer and 



Table 4. Comparison of the tidal radius of each dbh system at 
pericentre passage for the different components. Column (1) gives 
the model name, (2) production (Y) or not (N) of a tidal tail. Col- 
umn (3) gives the ratio between the tidal radius (r-pi) versus the 
radius enclosing 99% of the luminous mass (rperii) of the smallest 
system at pericentre. Column (4) lists the same information for 
the largest system. The tidal radius is computed with the impul- 
sive approximation. 



model 


Tail 




7^T9 /?'nnri9 
' ± ^ / ' ^ UI 1 z 


(1) 


(2) 


(3) 


(4) 


ip 


Y 


0.912 


1.117 


la 


Y 


1.142 


1.107 


Ir 


N 


1.096 


1.134 




Y 


0.236 


0.292 


2p 


N 


2.010 


2.264 




Y? 


1.258 


1.695 




Y 


0.341 


0.804 


2a 


N 


2.096 


2.219 




N 


1.557 


1.700 




Y 


0.504 


0.865 


2r 


N 


2.105 


2.383 




N 


1.369 


1.733 




Y 


0.439 


1.075 


3p 


N 


1.990 


2.302 




Y? 


1.231 


1.757 




Y 


0.300 


0.936 


3a 


N 


1.980 


2.251 




N 


1.427 


1.809 




Y 


0.500 


1.057 


3r 


N 


2.020 


2.330 




N 


1.153 


1.693 




Y 


0.257 


0.997 



more sharply defined. But extended tidal material is also 
present in r models, in which both spins are retrograde with 
respect to the orbit. 

Clues to this behavior come from the analysis of which 
precursor galaxy hosted the stars now in each tidal tail. Such 
information is also useful for the interpretation of observed 
merger remnants. We find that all the tails in our non-equal 
mass merger remnants originate from the less massive sys- 
tem. In particular, this is true for the a models, in which 
the less massive galaxy's spin is antiparallel to the orbit, 
and the primary spin is aligned with the orbit (cf. Table|2|l. 
Clearly, the tails in the a systems do not arise from spin- 
orbit coupling. Rather, tails in our strongly-interpenetrating 
merger orbits arise from tidal impulses. Wherever the tidal 
radius is smaller than the radius of the disc, material must 
be expelled from the galaxy. We estimate the strength of the 
tidal impulse using the impulsive approximation (e.g. Bin- 
ney & Tremaine 1987) to compute tidal radii. Table 0| lists 
the ratios between tidal radius and disc radius, measured at 
each pericentre passage, for all the dbh models. We also list 
whether tails appear at each given pericentre passage. (In 
equal-mass encounters, long tails form in the first passage 
through the pericentre, with subsequent pericentre passages 
yielding additional, shorter tails that sometimes project on 
to the longer ones; in non-equal mass models, tails form in 
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Figure 5. Time evolution of the luminous matter in run Ip of the disc-bulge-halo series. Time is indicated in the upper right corner. 
The half mass radius of the final merged system is 0.52. 



the second or third pass though pericentre). Table |1] shows 
that tails appear whenever rx/rperi < 1, showing that the 
tidal impulse is indeed responsible for the formation of tails. 
This explains that tails appear in models lacking spin-orbit 
coupling. Only in equal-mass mergers does the recipe fail, 
possibly because the impulsive approximation fails in this 
case. Tidal impulse tails are also formed from the large sys- 
tem, for example in models 3p and 3r. This tails are weak 
and do not extend much. 

The fraction of mass expelled into tidal tails varies with 
the mass ratio and the spin-orbit coupling. It ranges from a 
mere 2% for model 3r to up to 10% for Ip . Of this, a large 
fraction has negative binding energy, and is falling, or will 
fall back, on to the remnant (Hibbard & Mihos 1995). Only 
in model Ir do we find the tail material with positive energy 
(all particles in the tail at radius larger than 20 units: about 
1% of the luminous mass of the total system). This small 
escape fraction is a result of the sub-parabolic nature of our 
merger orbits; the escape fraction is 6% in the models of 
Barnes (1992). 



Rather different results hold for dh (bulge-less) systems. 
Figure |H1 shows snapshots of the large-scale distribution of 
luminous matter at the end of the nine dh simulations, as 
viewed from the z-axis. Tails appear only in models 2p* 
and 3p* (prograde, unequal mass mergers). These tails are 
composed of secondary material, indicating that spin-orbit 
coupling is the dominant mechanism for tail formation in 
these models. In general though, comparison with Figure |3 
reveals a general lack of extended tidal material in db merg- 
ers as compared with dbh models; this is surprising if we 
recall that mass ratios, merger orbits and spin orientations 
are identical in both sets of models. By inspection of the 
luminous mass distributions at intermediate stages during 
the merger orbit, we find that weak tails did actually de- 
velop in most of the models; however, this tail material has 
already fallen back to the main body by the time the merger 
is complete. 

The strength of tidal features therefore depends not 
only on tidal strength and spin-orbit coupling but also on 
the presence of a bulge in the precursor galaxies. This is 
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due to bar formation in the merging discs. As highlighted 
by Hernquist (1992), the central bulge in dbh models helps 
stabilize the disc against bar distortion. Bulge-less models 
develop strong bars due to the perturber's tidal field. The 
interaction of these bars with the haloes drives angular mo- 
mentum from the discs to the haloes, depriving the discs 
from the angular momentum needed for the development of 
large tails. Bar formation occurs in the bulge-less models of 
Hernquist (1992); however, some of his initial configurations 
have a perfect alignment of the discs with the orbital plane, 
creating a maximal spin-orbit coupling; this might lead to 
the formation of more prominent tails than we find here. 

The strong disparity in tail strength, geometry and life- 
times between our dbh and db models, (cf. Figs. |7| and (HJ 
highlights the difficulties in constraining the masses and po- 
tential depths of the dark matter haloes of galaxies from the 
lengths of the tidal tails, as attempted by Dubinski et al. 
(1996), Mihos et al. (1998) and Dubinski et al. (1999). In this 
last work the authors perform a survey to check the infiuence 
of halo mass profiles on the development of tails. They give 
the quantitative relation iic/wc 2.5 to produce tidal tails. 



with — 2|(/!>(r)| the escape velocity and v1 = GM{r)/r 
the circular velocity at a radius of approximately two times 
Rd- Our initial models do show different concentrations for 
the halo potential so we could expect differences in tail pro- 
duction. For our initial disk models the values for the Va/vc 
ratio are 2.12 for the dbh model and 2.19 for the dh case, 
which are both in agreement with the previous relation, and 
therefore we would not expect the differences reported ear- 
lier. Still, different tidal tails are produced, thus highlighting 
the role played by the bulge in stabilizing the disk against 
bar formation. 

3.2 Shells 

Hernquist & Spergel (1992), first noted that major mergers 
of disc galaxies may lead to the formation of shells, and we 
find that this is the case in our models. We extend their 
results by analysing mass ratios different than unity and 
model galaxies with and without bulges. On the other hand, 
our lower number of particles results in a lower definition for 
outlining existing shells. 
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Figure 7. Final distribution of the luminous matter in eacli of the dbh simulations. Each box measures 1.4 Mpc on a side, or approx. 
200 effective radii Rco{ the final models. The frame sizes are chosen to comprise out to the outermost tidal material. The systems are 
viewed from the positive 2:-axis; the merger orbits run counter-clockwise on the plane of the figure. 



The presence of shells in our merger remnants may be 
seen in Figures l^and lTHl which provide close-up views of the 
luminous mass distribution for dbh and dh merger models, 
respectively. These Figures confirm that shells do form in 
S-S mergers, without recourse to the technique employed by 
Hernquist & Spergel of "filtering" the particle distribution 
to plot only particles with near-zero radial velocity, i.e. near 
their apocentres. From the particle distributions in Figures|5] 
and IIUI the numbers and sharpness of shells vary with the 
mass ratio and with the presence or absence of a central 
bulge. Shells are much less prominent in dbh models (Fig.^Il, 
where we find shell-like structures in models Ir, 2r and 3a. 
Sharp features are present in models Ip and la, but those 
link with tails further out and may be described as well as 
inner, wound-up tails. 

Models without bulges lead to more prominent shells 
(Fig. llUH . Shells, of type II (all-round) or type III (irregular) 
following the definitions by Prieur (1990), are seen in models 
2a*, 2r*, 3a*, 3r*. We note the lack of shells in remnants 
of equal-mass mergers and on all prograde mergers. This 
contrasts with the shell system presented by Hernquist & 



Spergel (1992), a prograde merger of two equal-mass, bulge- 
less discs. The perfect alignment of the disc spins with the 
orbital angular momentum may have favored the formation 
of shells in their model. 

The lack of particle resolution prevents further anal- 
ysis of the shell distribution, radial range or interleaving. 
We can nevertheless enquire whether the shells are phase- 
wrapped or space-wrapped, by switching from coordinate 
space to phase space. Figure [TTI shows four representations 
of phase space for models 2a and 2a*. In Vi - r space, the 
V-shaped features provide the classical signature of phase- 
wrapped shells from material on radial orbits (Quinn 1984; 
Merrifield & Kuijken 1998). Model 2a shows very faint (if 
any) shells, while these are more prominent for 2a*. Whether 
these shells are indeed on radial orbits may be investigated 
by plotting K and Vt vs. the energy E (Fig. 1111 right pan- 
els). In Vi vs. E, vertical bands trace phase- wrapping, i.e. 
material with a range of radial velocities but common en- 
ergy and apocentre distance. For model 2a*, at the energies 
of the three most prominent bands, the Vt vs. E diagram 
shows clumps near K = 0. Therefore, the strongest shells 
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Figure 8. Final distribution of the luminous matter in each of the bulge-less simulations. Each box measures 1.4 Mpc on a side, or 
approx. 200 effective radii i?oof the final models; this is the same spatial scale as used in Figure |2l The systems are viewed from the 
positive z-axis; the merger orbits run counter-clockwise on the plane of the figure. 



in this model are radial, phase- wrapped shells. On the other 
hand, a clump of particles with K ~ 0.5 and = is 
seen (E « —0.4). This configuration may correspond to a 
space-wrapped shell. 



Sharp features in phase-space distributions are less 
prominent for the dbh models, thus confirming the seem- 
ingly inabiltiy of galaxies with prominent bulges to form 
shells through mergers. One could argue that the material 
in the tails could produce more prominent shells for dbh 
models when that material falls back to the main body. We 
have followed model 2a futher in time and found that there 
is no major difference with the picture given here, confirming 
the result given above. 



Gonzalez-Garcfa & van Albada (2005a, 2005b) find that 
major mergers between elliptical galaxies also produce faint 
shells when there is a contrast in the potential wells of the 
progenitors. 



3.3 Axis ratios 

We study the intrinsic shapes of the luminous mass distri- 
butions by computing axis ratios b/a, c/a, as was done in 
de Zeeuw & Franx (1991)'s study of the intrinsic shapes of 
elliptical galaxies. The quantities a > b > c are the axes of 
the inertia ellipsoid. For a homogeneous ellipsoid with axes 
2a, 26 and 2c, the eigenvalues of the inertia tensor are: 

El = {b^ + c^)/5, E2 = {a'' + c^)/5, E3 = {a +b'') /5. (10) 

Defining the axes such that E\ < E2 < E3,, it follows that 
the axis ratios are: 



E3 + E2 — El 
-Bi + -E3 — -E'2 ' 

E3 + E2 — El 
E1+E2-E3' 



(11) 



(12) 



We compute these ratios for radii including 30%, 50% and 
80% of the luminous mass. 

The axis ratios for dbh and dh models are plotted in 
Figure [T^ top and bottom, respectively. In such diagrams. 
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Figure 9. Close-up views of the final distributions of luminous matter for the nine dbh simulations. The size of each box is '~40 _Rc(see 
Table 151 . The systems are viewed from the positive 2-axis. 



spherical objects {a = b = c) lie at the upper right corner, 
oblate objects lie on the b/a = 1 line, and prolate objects 
lie on the diagonal {b/a = c/a). 

Figure fT?l shows a strong dichotomy of shapes between 
dbh and dh models. Mergers with bulges yield remnants 
that are oblate at all radii, whereas dh models are strongly 
triaxial, approaching prolate shapes, out to the inner half- 
mass radii. Further out, dh models approach the oblate line, 
a trend toward outer oblate shape already noted by Gerhard 
(1983), Barnes (1992) and Hernquist (1992, 1993) for equal- 
mass mergers. 

The triaxiality of dh merger remnants may be under- 
stood from the susceptibility of bulge-less galaxies to become 
bar-unstable during orbital decay, due to the tidal field: by 
the time the two cores are about to merge, their shapes 
are nearly prolate. In contrast, thanks to the stabilizing ef- 
fect of the bulge potential, dbh models are robust against 
bar formation, and the discs as well as the bulges remain 
nearly oblate throughout orbital decay. Thus, bar instability 
of bulge-less discs, which explained tail morphology differ- 
ences for dh and dbh merger remnants ('ii l3.lL lies also at 



the origin of the differences in final shapes between dh and 
dbh models. 

Franx et al. (1991), and Statler (1995) showed that real 
ellipticals have nearly oblate shapes. In the light of Fig- 
ure 1121 such intrinsic shapes are hard to explain within 
the coUisionless merger hypothesis for ellipticals, unless the 
progenitor discs had prominent central bulges before the 
merger: bulge-less galaxies generate remnants that are too 
triaxial. 

For dbh models, the most significant trend is an outward 
increase of c/a fiattening (e.g. from c/a — 0.7 at 30% radius 
to c/a — 0.55 at 80% radius, for model Ip). This is most 
likely due to the inner axis ratios reflecting the shapes of the 
merged spheroids (the bulges), while the outer distributions 
reflect the scattering of disc material over the merger orbital 
plane. Equal-mass dbh models also show mild triax;iality in 
their inner regions {Ip, Ir, la, 30% mass diagram). This may 
be explained by the general rule that equal-mass galaxies 
embedded in extended haloes merge on nearly-radial orbits, 
due to absorption of orbital angular momentum by the halo 
(e.g. Barnes 1992, Gonzalez-Garci'a & van Albada 2005b). In 
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Figure 10. Close-up views of the final distributions of luminous matter for the nine dh simulations. The size of each box is ~40 _Ro(see 
Table 13 . The systems are viewed from the positive z-axis. 



unequal mass mergers, in contrast, the smaller galaxy spirals 
in around the dominant galaxy, and the final shape is more 
nearly oblate. 

We are unable to reproduce Barnes' (1999) claim that 
3:1 mergers consistently lead to more fiattened ellipticals 
than equal-mass mergers. Naab & Burkert (2003) have 
shown that 3:1 and 4:1 remnants have higher ellipticities 
than 1:1 remnants (see their figure 4). Figure 1121 shows 
that final flattening c/a does not scale with the mass ra- 
tio of the merging galaxies. And, in the cores (30% enclosed 
mass), equal-mass mergers actually lead to more fiattened 
structures, as a result of the head-on final merger discussed 
above to explain the mild triaxiality of these cores. Barnes' 
and our results may difi'er due to differences in the surveyed 
merger orbits, in the relative orientations of the initial spins, 
and in the bulge-to-disc ratios {B/D = 1 : 3 in Barnes', 
B/D = 1 : 2 in our models). We get similar axis ratios 
for 1:1 mergers (e.g., models 2-3-5 in Barnes 1992), argu- 
ing against systematic differences in the derivation of axis 
ratios. Both Naab & Burkert and Barnes have simulated 
larger samples than presented here for less massive bulges. 



Therefore the initial orientation of the disks is unlikely to 
result in the differences reported here. In either event, our 
models show that the flattening of 3:1 merger remnants is 
similar to that of the 1:1 remnants; this occurs whether or 
not the parent galaxies harbor central bulges. As a result, 
stating from axis ratio considerations that 3:1 mergers of 
disc galaxies lead to flattened, SO-type galaxies is not sup- 
ported by our data; the statement might be valid for higher 
mass ratios, 5:1 or 10:1. 



The above analysis indicates that flgure shapes and axis 
ratios of merger remnants strongly depend on the presence 
of a central spheroidal bulge in the precursor galaxies. Disc 
galaxies with central bulges yield oblate remnants similar to 
real ellipticals, while galaxies without bulges yield objects 
that are triaxial within the effective radius. That these ob- 
jects are rare among elliptical galaxies may be a new prob- 
lem for the merger hypothesis for the formation of ellipticals 
from collisionless mergers of bulge-less disc galaxies. 
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Figure 11. Left panels: Vi vs. R diagram for (top) simulation 2a, and (bottom) simulation 2a*. Right panels: radial and tangential 
velocities vs. energy for (top) simulation 2a and (bottom) simulation 2a*. 



3.4 Surface density profiles 

We derive surface density profiles from ellipse fits to the iso- 
density contours of the merger remnants, as viewed from a 
point of view parallel to the initial orbital angular momen- 
tum. Surface density profiles for the nine dbh models are 
shown in Figure ITsl while those of the dh models are given 
in Figure ITU 



These profiles are fitted using a least-square method to 
a de Vaucoleurs (1958) R^^^ profile: 

/ = Je exp {~7.61[{R/Rcf^'^ - 1]}, (13) 
and also to a Sersic (1968) i?^/" profile: 
/ = Je exp {(-1.999271 + 0.3271)[(i?/i?o)^/" - 1]}. (14) 
where R^ is the radius enclosing half the luminous mass, 
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Figure 12. Axis ratios of the final models, measured within radii 
enclosing (left) 30%, (middle) 50% and (right) 80% of the lumi- 
nous matter. Top panels are dbh models, bottom panels are dh 
models. The diamond indicates model 2aG. 

Ic is the luminous surface density at Re, and n is the profile 
shape index, with n — 1 corresponding to an exponential 
profile, and n = 4 to the de Vaucouleurs law. In order to 
avoid softening effects at small radii, we have only fitted 
points outside 2e, with e = 0.02 for dbh models and 0.09 for 
dh models. Fits extend out to 7? = 5, the truncation radius 
of the initial discs. Best-fitting parameters for the Sersic and 
de Vaucouleurs fits are given in Table |^ together with the 
rms from each fit. In Figures [T^ and ITU the bottom panels 
for each of the plots give the residuals from both fits. 

The surface density profiles are well approximated by a 
de Vaucouleurs profile, as found in previous works (Barnes 
1992, Hernquist 1992, 1993). However, inspection of the rms 
values (Table |SJ shows that a better fit is provided with the 
Sersic model, especially for the db models. This is reassuring, 
since the Sersic model provides a better fit to real ellipticals 
than the 7?^'''' profile as well (Caon et al. 1993). Remnants 
of dbh models show Sersic indices n = 3.3 to n = 8, while 
dh models have lower n values, n = 2.4 to n = 3.2. The 
lower n of the dh models are a result of the lower central 
densities of the initial models, highlighting that mergers of 
bulge-less discs cannot account for the central densities of 
elliptical galaxies (Hernquist 1992). However, recent results 
from Trujillo et al. (2004) show that giant elliptical galaxies, 
brighter than — 18A/s, cover a range of values for the Sersic 
exponent from 2 to 9, which would be in agreement with the 
range for dh experiments. We note that the Sersic indices of 
dbh models do not reflect the density profiles of the initial 
bulges, which are modeled as King profiles and therefore 
approach an exponential rather than a R^^* profile. 

All the profiles deviate from both the de Vaucouleurs 
and the Sersic laws at radii R > 3. The mass involved is 
small (~5% of the initial luminous mass). Most likely, the 
deviations reflect the initial disc truncation radius {R = 5) 



combined with the fact that the outer parts are not yet fully 
virialized: after merging, the systems were let to relax for 
about 10 time units (i.e. 40 half-light radius crossing times), 
which is of the order of the remnants' crossing times at a 
radius R = 3. 

Perhaps surprisingly, the Sersic index n does not de- 
pend on the mass ratio of the merged galaxies. Instead, the 
final profile shapes are affected by the alignments of the 
galaxy spins with the orbit. Prograde and anti-parallel con- 
figurations give higher values of n than orbits where both 
spins are opposite to the orbital angular momentum (see 
Fig. 1151 . This is true for both dbh and dh models, and proba- 
bly reflects more extended mass distributions resulting from 
higher absorption of orbital energy by the galaxy envelopes 
in prograde mergers. 

Because A''-body systems are dimensionless, the models 
cover a small base-line in luminosity (or mass) and are un- 
able to generate a luminosity- Sersic index n diagram that 
reproduces the trend followed by ellipticals (e.g. Caon et 
al. 1993). However, from the lower Sersic indices of bulge- 
less merger remnants, it appears that the Sersic index of 
the remnant is a tracer of the concentration of the precur- 
sors. We use here the Concentration Index Cro introduced 
by Trujillo, Graham, & Caon (2001). Cro is defined as the 
ratio of the fight within O.SRcto the Hght within Rc, where 
f?eis the effective radius derived from the Sersic fit. Crc is a 
monotonic function of n (Fig. 1151 top-right panel). The con- 
centration index Crc is distinctly higher for dbh models than 
for dh models. Clearly, the presence of the bulge in the pro- 
genitors leads to a higher central concentration in the light 
distribution of the remnants. Higher central concentrations 
are associated to higher central velocity dispersions (Fig. 1151 
bottom-right panel), as found in real ellipticals (Graham, 
Trujillo, & Caon 2001). 

3.5 Boxiness-diskiness 

We calculate the boxiness-diskiness of the isodensity con- 
tours of the merger remnants by fitting the density residuals 
5{4>) along the best-fitting ellipses to a Fourier series (Carter 
1978), 

S{(p) = 5 + ^ anCoa{n(j)) + ^ b„sin{n(l>), (15) 

where is the azimuthal angle sweeping the ellipse, 5((/!>) is 
the residual intensity along the best-fitting ellipse, S is the 
mean residual intensity and n = 1,2,.... For well behaved 
isophotes, S, an and b„ should be small. Isophotes are called 
disky when 04 is positive, and boxy when 04 is negative (e.g. 
Bender 1988). 

We show profiles of the 04/0 parameter for eight of our 
models in Figure [T^ The view point is parallel to the inter- 
mediate axis of the inertia ellipsoid, and the average of 60 
snapshots, closely spaced in time, has been taken to increase 
the signal-to-noise. Amplitudes of the boxiness/diskiness de- 
viations are a few percent, but show peaks reaching 10%, 
which is of concern as 04 deviations in ellipticals are typi- 
cally < 3% (Bender 1988). However, a number of patterns 
are apparent in Figure lTSl which suggest that we are not just 
measuring noise. We find that the isophotal properties of 
dbh and dh remnants are quite different, dbh models gradu- 
ally change from boxy (Ir) to disky (3p) as the mass ratio of 
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Figure 13. Surface density profiles for tlie nine dbli models. For each model, the top panel shows the surface density profile (dots), the 
best-fitting R^^^ model (dotted line), and the best-fitting Sersic model (dashed line). The bottom panel shows the residuals from the 
two fits. 
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Figure 14. Surface density profiles for the nine dh models. For each model, the top panel shows the surface density profile (dots), the 
best-fitting R^/'^ model (dotted line), and the best-fitting Sersic model (dashed line). The bottom panel shows the residuals from the 
two fits. 
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Table 5. Best-fitting parameters for Sersic and de Vaucouleurs fits to the surface density profiles. For each 
model, the first line gives the best-fit parameters for a Sersic profile, and the second line the parameters for 
a de Vaucouleurs fit (n = 4). 
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Figure 15. Top left: mass ratio versus index n of Sersic profile 
fit. Top right: Trujillo et al (2001) concentration parameter C'n^ 
versus n. Bottom left: Cjje versus mass ratio. Bottom right: Cn^ 
versus central velocity dispersion. Open symbols are dbh models 
while filled symbols are dh models. Circles are 'p' models, stars 
are 'a' models and triangles are 'r' models. The diamond indicates 
model 2aG. 



Models without bulges show an altogether new and dif- 
ferent pattern, as their main bodies are strongly boxy. Such 
isophotal profiles are not found in elliptical galaxies, and, to- 
gether with their strong triaxiality (5 13.311 represent another 
failure of coUisionless, bulge-less galaxy mergers to yield ob- 
jects resembling real ellipticals. 

Boxiness-diskiness has been claimed to be point-of-view 
dependent (Fasano 1991, StiaveUi et al. 1991, Ryden 1992, 
Governato et al. 1993, Heyl et al. 1994, Gonzalez-Garcfa 
2003, Naab & Burkert 2003). We have studied the point-of- 
view dependence of the 04 parameter of model 3r by plot- 
ting histograms of the mean values of 04 inside 1.5Rc and 
3Rc for one hundred random viewpoints (Fig. I17II . A time 
average over 60 snapshots was taken to improve the signal- 
to-noise ratio. Overall, the 04 parameter is found to vary 
from positive to negative, in agreement with previous works, 
although the distribution is often skewed (Fig. 1171 left panel) 
indicating a larger probability to show a distinct isophotal 
signature (boxy in the plotted example). Clearly, the distri- 
bution of 04 depends on the radial range of the measurement 
(cf. left vs. right panels in Fig. I17II . highlighting that com- 
parison to observations needs to be performed over similar 
radial ranges. 



the merging galaxies increases. Naab et al. (1999) reported 
the same trend, which prompted these authors to conjecture 
that intermediate-luminosity ellipticals might be remnants 
of mergers of unequal disc galaxies. Naab & Burkert (2003) 
also conclude that a trend from boxy to disky can be ob- 
served with the mass ratio. Whether such association is real 
is currently unclear; Naab & Burkert (2001) rejected it on 
the basis of the line-of-sight velocity distribution properties 
of the remnants. 



3.6 Rotation curves 

Whether mergers produce rotation curves with the ampli- 
tudes observed in elliptical galaxies is a key test on the 
merger hypothesis. We show major and minor axis rota- 
tion curves for all the merger remnants in Figure [T51 Rota- 
tion curves were derived by looking at the models from the 
2/- axis and using a viewing slit of length of 2 units and a 
width of 0.3 units; such rotation curves extend up to 2Rc in 
all cases. All the remnants show rotation along the major 
axis and nearly zero minor-axis rotation, i.e. they rotate like 



© 0000 RAS, MNRAS 000, 000-000 



Elliptical galaxies from mergers of discs 17 





1 r- 

1 r 




1 r 



1 2 



2a : 



2 3 4 5 

R/R, 



2 3 4 5 

R/R, 



3p : 



2 3 4 5 

R/R, 



2 3 4 

R/R, 



OX. 



R/R, 



R/R, 



Figure 16. Boxiness-diskiness, as measured by lOOai/a, vs. 
R/Rb for several models. The point of view is parallel to the 
intermediate axis at the half-light radius of the remnant. 




Figure 17. Histograms of boxiness distribution for model 3r. 
We have calculated the mean 04 parameter inside 1.5-Re(left) and 
3i?e (right) of this model as measured from one hundred points of 
view (randomly chosen) for 60 snapshots . 



oblate ellipsoids. In particular, this is true for most of our 
1:1 mergers, which are commonly assumed to lead to kine- 
matic misalignment (e.g. Heyl et al. 1996, Naab & Burkert 
2003). Of the 1:1 mergers, only model la shows a mild mis- 
alignment. A slight nuclear minor axis rotation is also seen 
in model 3r. 

Two factors affect the rotation amplitude. First is the 



spin alignment with the orbit: as expected, p models show 
the highest rotation. Second is the presence of a bulge in 
the progenitors: dbh remnants show both a higher rotation 
amplitude than their corresponding dh models, and a dis- 
tinct rotation profile shape, given by a rapid inner rise and 
an outer plateau. In contrast, the rotation curves of dh rem- 
nants show a gentle, uniform rise to the outermost mea- 
sured point. Clearly, the bulge components, which remain 
self-gravitating during most of the merger, deliver their or- 
bital angular momentum further in than bulge-less galaxies. 
These results are summarized in Figure lT^ which plots Vmax 
for all the models vs. the mass ratio. Model 3r* is the only 
model showing no rotation; this 3:1 retrograde merger of two 
dh galaxies shows nearly zero rotation, on both the major 
and minor axes, and a mild nuclear counter-rotation. Its con- 
figuration is similar to those studied by Balcells & Gonzalez 
(1998) to produce counterrotating cores from mergers of disc 
galaxies. 

We obtain the degree of rotational support by normaliz- 
ing the rotation curves by the velocity dispersion. In Figure 
1201 we plot V{R)/a{R) vs. R/Rc for each model, where R is 
the projected radius. As before, these plots are produced for 
a point of view parallel to the y-axis. A mean has been done 
for both sides of the rotation curve. The patterns outlined 
above for V{R) are reproduced here, with V{R)/a{R) being 
higher for dbh remnants and for direct configurations. Addi- 
tionally, Figure I2S1 shows that the amplitude of V{R)/<j{R) 
is higher for higher mass ratios (the dependency of the ro- 
tation with mass ratio could not be derived from Figure [TSl 
as the total mass varies, from 2 to 4, between 1:1 and 3:1 
models). This effect mostly results from the lower velocity 
dispersions in the outer parts of 2:1 and 3:1 remnants as 
compared to the 1:1 remnants, i.e. lower-mass companions 
produce less heating in the outer parts of the larger disc, 
and its cold dynamics is partly retained in the remnant. 

The variety of rotation amplitudes, for given mass ra- 
tio and spin orientation, emphasizes that the rotation of 
merger remnants arises from a combination of retained ini- 
tial spin and redistribution of orbital angular momentum 
(Bendo & Barnes 2000). Therefore, ellipticals of merger ori- 
gin should show a wide range of rotation amplitudes. Fig- 
ure I2UI may be directly compared to observations. While 
giant ellipticals have long been known to rotate slowly, 
with Vmax/cFQ « 1 (lUiugworth 1977), recent, deep spec- 
troscopy of intermediate-luminosity ellipticals shows high 
rotation, with Knax/o" in the range of 1-3 (Rix et al. 1999; 
Vinax/o" ~ 3.8 for an SO galaxy in their sample). These val- 
ues have been compared to A'^-body 3:1 merger models by 
Bendo & Barnes (2000) and by Cretton et al. (2001), with 
conflicting conclusions. At i? = 2i?e, our 3:1 mergers yield 
at most Knax/cr = 2 (model 3p), hence we concur with Cret- 
ton et al. that 3:1 mergers do not yield the highest Knax/o" 
values observed in intermediate-luminosity ellipticals. We 
note however that, out to R ~-Re, our model 3p matches 
the fastest-rotating ellipticals. The discrepancy sets in at 
larger radii: although V{R)/a{R) continues to rise to the 
limit of our measurement, its amplitude at 1 < R/Rc< 2 
is lower than observations. The rotation curve y' G Mr/r 
is halo-dominated at these radii; more massive, extended 
haloes might yield faster rotation in the merger remnants. 
Alternatively, the coUisional formation of an outer gaseous 
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Figure 18. The left panels show the rotation curves of the luminous matter for the remnants of the 18 simulations, as measured by 
placing a viewing slit of width of 0.3 length units along the major axis. Models are viewed from the positive y-axis. The right panels 
show the rotation curves along the minor axis from the same point of view. The vertical bars indicate representative velocity error bars, 
computed using Poisson statistics. 



0.8 




1 1 1 

p — models 


1 1 1 


1 1 1 1 








— models 












r— models 








0.6 










^ 3p ; 


0.4 










..3a*- 


0.2 


- .1.0*- 








~---.3r* ' 


0.0 


, , 





, 1 , 


, 1 , , 


, , 1 , 


1.0 


1.5 


2.0 


2.5 


3.0 



Moss ratio 

Figure 19. Maximum rotation velocity, Vmax, vs. merger mass 
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disc, with ensuing star formation, could lie at the origin of 
the high outer Knax/o" values observed by Rix et al. 

3.7 Rotation support vs pressure support 

Whether eUipsoids owe their flattening to rotation or to 
anisotropy of the velocity dispersion tensor can be stud- 
ied with the Knax/o"o vs e diagram (Fig. . Oblate ellip- 
soids flattened by rotation show a well defined dependency 
of Knax/o"o with the ellipticity, given by the solid line in 
Figure 170 (Binnev 1978). In this diagram, giant elliptical 
galaxies fall well below the oblate line (lUingworth, 1977), 
while low-luminosity ellipticals and bulges scatter along the 
oblate line (Davies et al. 1983; Kormendy & lUingworth 
1982, 1983). The location of our merger remnants in the 
Knax/o"o vs £ diagram has been determined for 100 random 
view-points for each model. The results are shown in Fig- 
ure]^ (left panels), while the mean position for each model 
is plotted in Figure ITTI fright panels). 
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Figure 20. V {K) j a{K) vs. radius R/R^. Notation as in Figure 
1181 Top panels: dbh mergers. Bottom panels: dh mergers. Error 
bars for one of the highest R/Rcare given, computed using Pois- 
son statistics. 



The mean for bulge dominated models lies close to the 
isotropic oblate rotator line, dbh models do keep a high 
amount of rotation while being fairly round (e = 0.2). The 
spin configuration plays a role, as direct orbits yield higher 
rotation support than antiparallel or retrograde orbits. The 
mass ratio also plays a role since equal mass mergers are 
less rotationally supported than non-equal mass mergers, as 
noted by Naab et al. (1999) and Naab & Burkert (2003). 
Overall, the ranges in V/a and in ellipticity are similar to 
those of intermediate-luminosity ellipticals. 

Remnants of bulge-less systems show a larger scatter in 
the Knax/cro vs e diagram. Remnants of unequal mass galax- 
ies lie close to the oblate line, but equal-mass mergers lead to 
highly-flattened, slowly rotating systems that are pressure- 
supported. As discussed in 5 13.31 these systems are strongly 
triaxial as a result of transient bars during the merger evo- 
lution. The overall distribution in the Vmax/co vs e diagram 
is different both from that of intermediate-luminosity ellip- 
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Figure 21. Model rotation, expressed as Vm/cro vs. the projected 
ellipticity. Top panels: dhh mergers. Bottom panels: dh mergers. 
Measurements are made for one liundred points of view for eacli 
model, and are plotted with small dots in the left panels. The 
mean for each model is plotted with the model name in the right 
panels. Observations from Davies et al. (1983) are shown in the 
left panels; open circles are bright ellipticals, filled circles are low 
luminosity ellipticals and crosses are bulges. 

ticals and that of giant ellipticals, providing another line 
of evidence that mergers of bulge-less galaxies do not yield 
objects resembling elliptical galaxies. 

3.8 Anisotropy 

We have calculated the anisotropy parameter /3 as a function 
of radius for our models. This quantity is defined as: 

(16) 

where ai is the radial and at the tangential velocity disper- 
sion: 

= aj + al. (17) 

Radial anisotropy yields /? > 0, while /? < corresponds 
to tangential anisotropy. Profiles of /3 vs. radius are shown 
in Figure 1^ The top panel shows the dbh models, while dh 
models are shown in the bottom panel. Clearly, most rem- 
nants show tangential anisotropy. This is most pronounced 
in dbh remnants, which have tangential anisotropy at all 
radii. Anisotropy in strongest for unequal-mass remnants, 
an indication that these systems retain a dynamical mem- 
ory of the disc kinematics of the larger precursor disc. The 
only models deviating from tangential anisotropy are model 
la, which is nearly isotropic over the entire radial range, and 
model Ir which is radially anisotropic. This is likely a re- 
sult of the cancellation of spin and orbital angular momenta 
which leads to dominantly radial motions in the remnant. 

Remnants of dh models show more diversity in their 
anisotropy characteristics. All models show some degree of 
radial anisotropy inside the half mass radius. Such pattern is 
related to bar formation during the merging stages, as men- 
tioned in S 13.31 In the outer parts (r/ri/2 > 1), similarly to 
dbh remnants, most dh systems are tangentially anisotropic; 
some are isotropic, while spin cancellation leads model Ir* 
to show radial anisotropy. 
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Figure 22. Anisotropy profiles (3 vs. radius. Top panel: dbh mod- 
els. Bottom panel: dh models, 'p' models are depicted with solid 
lines, 'a' models with dotted lines, and 'r' models with dashed 
lines. The thickness of the lines gives the mass ratio; the thinnest 
lines are for 1:1 mass ratios and the thickest for 3:1. 

The tangential anisotropy characteristics are at odds 
with common lore that elliptical galaxies should be radially 
anisotropic at large radii (e.g., Berlin, Saglia, & Stiavelli 
1992, Berlin & Stiavelh 1993). This belief is rooted in gen- 
eral properties of violent relaxation. Under strictly violent 
relaxation conditions, stars lose their dynamical memory 
and owe their final dynamics to the energy exchanges with 
the time-varying global modes of the potential. Because in 
mergers 'all the action' occurs at the centre of the system, 
stars at large radii must have acquired their energy near 
the centre and henceforth must be on dominantly radial or- 
bits. The tangential anisotropy observed here is at odds with 
this framework, and emphasizes that merger relaxation is 
far from being completely violent (Zurek, Quinn, & Salmon 
1988). Owing to the preservation of the energy and angular 
momentum ordering during the merger, the outer parts of 
the remnants are comprised of stars in the outer parts of the 
discs, the tangential anisotropy being a relic of the disc kine- 
matics of the larger galaxy. Such preservation is maximal in 
dbh mergers, due to the stabilizing effect of the bulge's po- 
tential. In our models, radial anisotropy is restricted to two 
situations: when spins and orbital angular momentum cancel 
out in retrograde, equal-mass encounters, or in the central 
parts of dh remnants, as a result of transient bars triggered 
by the tidal field; in neither case is the radial anisotropy 
related to violent relaxation phenomena. 

Winsall & Freeman (1993) do not find evidence for 
anisotropy in giant ellipticals. Bender, Saglia & Gerhard 
(1994) study 44 elliptical galaxies fitting Gauss-Hermite 
polynomials to the LOSVD for each of them. The coeffi- 
cient of the four order polynomial is related to the tangen- 
tial anisotropy. Bender et al. do not find a clear trend in 
anisotropy, but find evidence that more rotation means more 
positive /i4 (negative values would mean radial anisotropy), 
although the conclusion is not firm. 

4 EFFECTS OF PARTICLE NUMBER 

As described in Section model 2a was replicated at high 
resolution (5 times more particles). Here we compare the 
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high-resolution model {2aG) and its low-resolution counter- 
part (2a) over the set of properties analysed in previous 
sections; resolution limitations could in principle affect all 
of them. The high-resolution model should be more accu- 
rate, therefore the differences between models 2a and 2aG 
provide a quantitative estimate of the errors in the models 
due to low particle number. 

Of particular concern was the ability to generate shells 
in \ow-N models due to discreteness effects (Hernquist & 
Spergel 1992). Like all dbh models, model 2a had been 
shown to have small, if any, shells (Figure O, which led 
to the claim that mergers of galaxies with central bulges do 
not form prominent shells (Section 13.211 . For run 2aG, sharp 
structures in the projected particle distribution are similar 
to those of model 2a. In phase space, where shells are gen- 
erally easier to identify, the distribution for model 2aG is 
quite similar to that of model 2a shown in Figure 1111 In 
particular, only the clumps at low binding energy that show 
up in Figure 1111 are visible in the corresponding diagram 
for model 2aG. This similarity suggests that shells are in- 
deed faint or non-existent for mergers of galaxies harboring 
central bulges. 

Variations in structural parameters are also small. Axis 
ratios are smaller by 9% at the 30% mass radius, and by less 
for outer radii; see Figure IT^ where model 2aG is displayed 
with a diamond. The lower mass ratios may reflect a stronger 
inward transport of orbital angular momentum as the nuclei 
of the high-resolution models manage to get closer to each 
other before merging together. The surface brightness profile 
of model 2aG is nearly identical to that of model 2a. The 
effective radius for model 2aG is 0.5% larger than that from 
run 2a. The Sersic exponent {n2aG = 4.82) is just 7% lower, 
while typical uncertainties for the Sersic exponent are of 
order 10%; and the concentration parameter is 2.5% lower; 
see Figure [TKl there model 2aG is shown as a diamond. 

The rotation profile from model 2aG is shown in Fig- 
ure l20l top middle panel, together with that from model 2a. 
The two curves are very similar. Differences arise in the in- 
ner parts, where particle resolution is more important. In 
the high resolution case, the V{R)/a{R) values are slightly 
higher at small radii, pointing perhaps to a more efficient 
transfer of the angular momentum to the inner parts. This 
may affect our results for other runs, pointing to an even 
steeper rise of the rotation curve for the inner parts, which 
would bring our models to a closer agreement with observa- 
tions (see Section f3.7ll . 

In conclusion, most properties of the simulations are un- 
changed with the increase of resolution, which suggests that 
particle discreteness effects are not blurring our principal 
conclusions. Some minor differences are observed, specially 
in the inner parts where particle resolution is better. 



5 DISCUSSION 

Our models show that the remnants of coUisionless mergers 
of disc galaxies have markedly different properties depend- 
ing on whether the precursors harbor central bulges. Mod- 
els with bulges (dbh) show higher central densities, oblate 
shapes, oblate-rotator kinematic support, and a mixture of 
disky and boxy isophotal shapes. Bulge-less models (dh) 
show shallow central densities, strongly prolate or triax- 



ial figure shapes, sub-oblate kinematic support, and boxy 
isophotes. In addition, dbh merger remnants show promi- 
nent tidal tails and no shells, while dh remnants do not have 
tails but show shells. 

We trace the origin of these differences to the fact that 
the tidal field distorts bulge-less discs into transient bars 
(Hernquist 1992). These remove angular momentum from 
the luminous bodies (Hernquist & Weinberg 1992), result- 
ing in low rotation, and generate prolate figures which, upon 
merging, yield prolate or triaxial remnants. The nearly ra- 
dial orbit makeup of tidally-distorted dh models may ex- 
plain the formation of phase- wrapped shells. In contrast, 
central bulges help stabilize the discs against bar distor- 
tion (Hernquist 1993), which allows the discs to retain part 
of their spin throughout the merger decay. This yields rem- 
nants with near oblate shapes and oblate-like kinematic sup- 
port (Barnes 1992, 1999; Bendo & Barnes 2000; Cretton et 
al. 2001). It also allows for more efficient spin-orbit coupling 
which is at the origin of tidal tail formation in prograde 
encounters (Toomre & Toomre 1972, Barnes 1992). 

Besides the well known difficulty in reproducing the cen- 
tral densities of elliptical galaxies (§0, dh merger models 
show additional properties which make them unlikely pre- 
cursors of today's ellipticals: they yield figure shapes that 
are too strongly triaxial, and their remnants show too lit- 
tle kinematic support. One could argue that their level of 
rotational support, their triaxial figure shapes and their 
boxy isophotes, make them akin to giant ellipticals; however, 
bulge-less galaxies are significantly less massive than galax- 
ies with massive bulges (e.g. Trujillo et al. 2001; Balcells, 
Graham, & Peletier 2004), and their luminosities are lower 
than those of giant ellipticals by factors of 10-100, hence an 
evolutionary link between dh mergers and anisotropic giant 
ellipticals cannot be established. 

Merger models with central bulges compare better with 
real ellipticals, specifically with intermediate-luminosity el- 
lipticals. In agreement with Hernquist (1993), our dbh mod- 
els yield realistic central densities. Further, figure shapes are 
oblate, and are supported by rotation. The level of rotational 
support, out to ~ i?e, compares well with even the fastest- 
rotating ellipticals. Note that, because matching the rotation 
requires disc kinematics in the precursors, and matching the 
central densities requires a bulge, precursors are required 
to harbor both a disc and a bulge if, as claimed by cur- 
rently popular galaxy formation models, ellipticals formed 
by mergers. 

A number of shortcomings are present in dbh coUi- 
sionless merger models. As noted above, the models can- 
not reproduce the rapid rotation in the outer parts of 
intermediate-luminosity ellipticals. Also, the way dbh mod- 
els reproduce the isophotal properties of ellipticals is not 
fully satisfactory. We concur with previous works (Naab et 
al. 1999; Bendo & Barnes 2000) that equal-mass mergers 
yield boxy isophotes, while 3:1 mergers yield diskiness. How- 
ever, not all initial orbit configurations of our 3:1 mergers 
produce disky isophotes; the values of the diskiness param- 
eter a4/a are higher than generally observed; and, a partic- 
ular system can be seen as either boxy or disky depending 
on the point of view. More generally, diskiness is a property 
of intermediate-luminosity ellipticals, and boxiness of giant 
ellipticals (Bender 1988; Rix & White 1990), a relation our 
models cannot explain given that gravitational N-body mod- 
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els are scale-free. The role of post-merger gas infall needs to 
be taken into account to explain the systematics of boxiness- 
diskiness with galaxy luminosity (khochfar & Burkert 2004) . 

On the whole, our models are more successful in re- 
producing oblate-rotator, intermediate-luminosity ellipticals 
than giant ellipticals, in general agreement with Naab & 
Burkert (2003). This may be surprising given that triaxial- 
ity and low rotation are properties generally expected from 
coUisionless merger dynamics. Pre-existing discs are a prob- 
lem for the formation of slowly-rotating, triaxial ellipticals, 
as the precursors retain too much of their spin and yield 
remnants that rotate too fast. Therefore, galaxies lacking or 
with less massive discs are potential merger precursors of 
giant ellipticals, as in the elliptical-elliptical merger models 
of Gonzalez-Garcia & van Albada (2005a,b). These may be 
seen as the final stages of the sequence of mergers galax- 
ies are believed to follow (Khochfar & Burkert 2003, van 
Dokkum & Ellis 2003, Bell et al. 2004). 
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